Broken translational symmetry in an emergent paramagnetic phase of graphene 
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We show that the spin-density wave state on the partially filled honeycomb and triangular lattices 
is preempted by a paramagnetic phase that breaks an emergent Z4 symmetry of the system, asso- 
ciated with the four inequivalent arrangements of spins in the quadrupled unit cell. Unlike other 
emergent paramagnetic phases in itinerant and localized-spin systems, this state preserves the rota- 
tional symmetry of the lattice but breaks its translational symmetry, giving rise to a super-lattice 
structure that can be detected by scanning tunneling microscopy. This emergent phase also has 
distinctive signatures in the magnetic spectrum that can be probed experimentally. 



I. INTRODUCTION 



Unconventional paramagnetic phases are characterized 
not only by the absence of long-range spin order, but also 
by a broken symmetry related to new degrees of free- 
dom that emerge from the collective magnetic behavior 
of the system. As a result, their elementary excitations 
and thermodynamic properties are rather different than 
those of an ordinary paramagnet. These phases usually 
appear in frustrated systems with localized spins, as a re- 
sult of the interplay between frustration and fluctuations. 
Canonical examples include the Ising-nematic phase of 
the extended Heisenberg model on the square lattice 
[]|, the spin-nematic phase of the Heisenberg model on 
the kagome lattice 0|, and the magnetic-charge ordered 
phase in kagome spin ice J3|. Itinerant magnetic sys- 
tems can also display paramagnetic phases with unusual 
broken symmetries. This is believed to be the case in 
the ruthenates [j| and in the iron-based superconduc- 
tors [5|-|7|. In these systems the emergent paramagnetic 
phase breaks the lattice rotational symmetry, while the 
spin-rotational and lattice translational symmetries re- 
main preserved. 

In this paper, we present an unusual itinerant para- 
magnetic phase that breaks the translational invariance 
without changing the point-group symmetry of the lat- 
tice. This phase arises in partially filled hexagonal (trian- 
gular and honeycomb) lattices, preempting a spin-density 
wave (SDW) order, and could potentially be realized in 
single-layer graphene doped near the saddle point of the 
band-structure (3/8 or 5/8 filhng) [ii|. The SDW order 
below T/v for fermions on a hexagonal lattice is uniaxial, 
with all spins pointing along the same direction |10| . The 
magnetic unit cell contains eight sites, six of which have 
moment —A and two have moment 3A, see Fig.[T] This 
state breaks not only the 0(3) spin-rotational symmetry, 
but also a discrete Z4 symmetry related to the four in- 
equivalent choices for the positions of the large 3 A spin 
moments in the eight-site unit cell. These four inequiv- 
alent spin configurations transform into each other upon 
translation of the origin of coordinates to neighboring 
hexagons - from point A to points B, C and D in Fig. [T] 




Figure 1: (Color online) The four inequivalent uniaxial SDW 
states. The quadrupled unit cell is highlighted with yellow 
boxes. Among the eight sites of the enlarged unit cell, two 
have large spin moment (Si) = 3A (red arrows) and six have 
small moment (Si) — —A (blue arrows). The total spin in 
each unit cell is zero. The four different states correspond to 
the additional Z4 symmetry of the order parameter manifold. 



Thus, breaking the Z^ symmetry corresponds to breaking 
the translational symmetry of the lattice. 

Of course, once the 0(3) symmetry is broken, the Z4 
symmetry has to be broken too. We show, however, that 
the Z4 symmetry breaks down at higher temperatures 
than the 0(3) symmetry. As a result, the SDW ordering 
at Tj\j is preempted by a phase transition at Tz^ > T/v, 
which falls into the universality class of the four-state 
Potts model. In the Z4 phase at Tn <T < Tz^, {S,) = 
for all sites (i.e., this phase is a paramagnet), and the 
unit cell is a hexagon (green dashed line in Fig. [3]), i.e., 
the Cg rotational symmetry of the lattice is preserved. 
Yet, the unit cell has eight inequivalent sites - for six of 




Figure 2: (Color online) Schematic representation of the emer- 
gent Z4 symmetry-breaking phase. The quadrupled unit cell 
is highlighted with yellow boxes; the numbers indicate the 
eight inequivalent sites. This state has the same broken trans- 
lational symmetry as that of the SDW in Fig. 1(a), but 0(3) 
symmetry is preserved (no arrows). The nearest-neighbor cor- 
relation function (Si • Sj) is A'^ for blue bonds and — 3A^ for 
red bonds. Other three states are obtained by moving the 
origin of coordinates from A to either B or C or D. 





Figure 3: (Color online) (a) The Fermi surface at the dop- 
ing level of interest is a hexagon inscribed within the hexag- 
onal Brillouin zone (BZ), for both honeycomb and triangu- 
lar lattices. The FS has three saddle points Ma located 
at the corners of the hexagon. Pairs of inequivalent saddle 
points are connected by three inequivalent nesting vectors 
Qi = (0,27r/y3), and Q2,3 = (=F7r/3, -tt/^). (b) Order- 
parameter space of the preemptive Z4 phase. In the ordered 
phase the vector order parameter = (i^i, (/>2, c/is) points to- 
ward one of the four corners of a regular tetrahedron. 



them (i = 1...6 in Fig. [2]) the bond correlators (S^ • Si+s) 
with their nearest neighbors are A^ and — 3A^ (blue and 
red bonds in Fig.[ljb)), while for the remaining two sites 
{i ~ 7,8) all bond correlators are — 3A^. The broken 
Z4 symmetry corresponds to choosing these two "spe- 
cial" sites out of the eight sites in the unit cell. One such 
choice is shown in Fig. [2l One can easily verify that the 
other three choices correspond to moving the origin of 
the coordinates from A to one of the points S, C, or D 
in Fig. [21 This obviously implies that the the transla- 
tional symmetry of the lattice is broken. Experimentally, 
the quadrupled unit cell in the Z4 phase can be readily 
probed by scanning tunneling microscopy (STM). Fur- 
thermore, we show that the transition to this phase is 
accompanied by a jump of the staggered spin susceptibil- 
ity, which can be probed by neutron scattering or nuclear 
magnetic resonance (NMR). 



II. THE UNIAXIAL SDW ORDER 

The Fermi surface (FS) of graphene near 3/8 or 5/8 fill- 
ing is near-nested and contains three saddle points with 
nearly vanishing Fermi velocity (the three Ma points in 
Fig.[2ta))- Pairs of inequivalent Ma points are connected 
by three commensurate nesting vectors Qi = (0, 27r/\/3) 
and Q2,3 ~ (=F7r/3, — 7r/-\/3). The divergent density of 
states at the M-points makes doped graphene a fertile 
ground for exploring nontrivial many-body density-wave 
and superconducting states (lll4l7| . The SDW instabil- 
ity is subleading to a chiral d-wave superconductivity ex- 
actly at 3/8 or 5/8 filling [ll|, but can become the leading 
instability slightly away from 3/8 or 5/8 filling |l(]|[l6| . 

In particular, the FS at the saddle-point doping lev- 
els, e.g. 3/8 or 5/8 for the honeycomb lattice, is a 
perfect hexagon inscribed within a hexagonal Brillouin 



zone (BZ) as shown in Fig. ^b). This FS is com- 
pletely nested by three wavevectors Qi = (0, 27r/-\/3), 
and Q2,3 — (±7r/3, — 7r/-\/3), and the nesting opens the 
door to an SDW instability. However, not all points on 
the Fermi surface are of equal importance. In particular, 
the three saddle points Ma (a = 1,2,3) give rise to a 
logarithmic singularity in the DOS and control the SDW 
instability at weak coupling. 

Thus, we consider the following Hamiltonian: 

H= ^ £aC'a^aCa,a (1) 

a=l,2,3 

where c^ ^ creates electrons with spin a around the sad- 
dle point Ma- There are two electron-electron interac- 
tions that contribute to the SDW channel, namely 32 
and 53, which represent the forward and umklapp scat- 
terings, respectively. The dispersions in the vicinity of 
the saddle points are 



3t 



ei(k) = ^{kl-3kl), 



3i 



£2,3(k) = --^2ky{kyTV3k,), 



(2) 
(3) 



where ii is the nearest-neighbor hopping constant. The 
quartic interaction terms in Eq. ^ can be decoupled via 
the Hubbard Stratonovich transformation with the SDW 



order parameters: A^ = Aa,i 
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Each of these vector order parameters corresponds to a 
nesting vector which connects two saddle points: Q.j = 
Ma — Mb. The partition function of the system can then 



be written as Z = J Vc^Vc DA, exp(-5[ct, c , A]) with 



-^T '-aj^a.a. 



S[c,c\A,]=Y,[cUd. 



(4) 
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where J^ = L dr. The fermionic part becomes 
quadratic and can be integrated out. By expanding the 
resulting action to fourth order in A^, we obtain the ef- 
fective action: 

5[A,]=ro5:|jA,p + ^|^(|Aip + |A,p + |A3p)^ 



(|Ai|2 + |A2p-2|A3p) +3(|Aip-|A2p) 



(Ai.A2)' + (A2-A3)2 + (A3-Ai)^ 



(5) 



Here tq oc (T — Tn), where T^ is the mean-field 
SDW transition temperature. The coefficients u,v,g 
in Eq. [5] were calculated in Ref. [10| and found to be 
positive, with v/u = f/log (PF/T/v) <C 1 and g/u = 
(Tn/W) I log (M^/Tjv) < 1, where W is the bandwidth. 

Minimizing S'[Ai] with respect to A^ and neglecting 
momentarily the fluctuations of the A^ fields, we see 
that u > implies that the magnitudes of A^ are equal, 
while 5 > makes all A^ collinear. The particular uni- 
axial state with (Ai,A2,A3) = (A,A,A)fi is shown 
in Fig. [BJa). There exists, however, three other states 
with the same energy, (A, —A, — A)n, (—A, A, — A)n, 
and (—A, —A, A)n. These states can not be obtained 
from the one shown in Fig. [IJa) by a global spin rota- 
tion. Instead, these four degenerate states are related 
by a translational Z4 symmetry - they transform into 
each other by moving the origin of coordinates from A to 
B, C, or D (Fig.[Ti;b)-(d)). The ground state in Fig.^a) 
chooses a particular direction of n and also one of the four 
positions of the origin of coordinates and therefore breaks 
0(3) X Z4 symmetry. 



III. PREEMPTIVE Z4 PHASE 

A. Order parameters 

We now allow A^ to fluctuate and analyze the possi- 
ble emergence of a phase in which Z4 symmetry is bro- 
ken but 0(3) symmetry is preserved. In such a phase 
(Ai) = 0, but (Ai • Aj) 7^ 0. A proper order parameter 
for the Z4 phase is the triplet (/> — (0i, 02,03), where 
^i = g {Aj ■ Ak) and (ijk) are cyclic permutations of 
(123). The Z4 symmetry breaking phase has (0i) = ±(f>, 
with the constraint (/)i0203 > 0. To investigate whether 
this state emerges we go beyond the mean-field approx- 
imation for S'[Ai] by including fluctuations of the A^ 



fields, and re-express the action in terms of the collec- 
tive variables 0^. We analyze this action assuming that 
fluctuations of 0^ are weak and check whether a non-zero 
((/)i) emerges above the SDW transition temperature. 

To obtain the action in terms of 0^ we apply a 
Hubbard-Stratonovich transformation [6| and introduce 
six auxiliary fields, one for each quartic term. These 
six fields include two fields ^1 ex (A^-I-Aj — 2A2) and 
C2 oc (A^ — A|) which break the Cq rotational symme- 
try, the three fields 0^ ex Aj • A^ associated with the Z4 
symmetry breaking, and the field ip ex (A^ + A2+A§) as- 
sociated with the Gaussian fluctuations of the A^ fields. 
Details of the Hubbard-Stratonovich transformation can 
be found in Appendix A. In particular, we show that the 
non-zero values of (^1 and C2 are energetically unfavor- 
able because w > so we set Ci = C2 = in the following 
analysis and consider states that preserves the lattice ro- 
tational symmetry. 

The quartic terms in Eq. ^ can be decoupled using 
the auxiliary fields (pi and ip. Because we allow the A^ 
fields to fluctuate, we include non-uniform space/time 
configurations, i.e., replace A — > A^.^ and rg — >■ rp -I- 
q^ + T\ujm\ hi Eq. ^, with Wm = 2'rmrT. Near a fi- 
nite temperature phase transition thermal fluctuations 
are the most relevant, and we restrict our analysis to the 
LUrn ~ component. The new action now depends only 
on the -0 and cp fields: 



S[^A] 
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2u 



log( 



dctA" 



(6) 



where |</,p = E.0f> /, = IJ^ I d'q, and V is the vol- 
ume of the system. The matrix X is 



X 




(7) 



with renormalized x 



To + ip + q'^ = r + q"^ 



absence of broken Z4 symmetry, long-range SDW order 
sets in at r = 0, hence an intermediate phase exists if Z4 
symmetry is broken at some r > 0. 

The action ^ is an unconstrained function of ijj, which 
is the usual situation for Gaussian fluctuations [20|, and 
reflects the fact that (Af ) ^ 0. However, we are princi- 
pally interested in the fields 0; , which have zero expecta- 
tion value in the absence of Z4 symmetry-breaking. The 
mean-field theory for the action (O is the set of coupled 
saddle-point equations - the minimum with respect to 
fluctuating fields 0^ and the maximum with respect to ijj. 



B. Mean-field theory 

The four possible realizations for the Z4 symmetry 
breaking correspond to 0^ — ±0 subject to the constraint 
0102 03 > 0. We substitute this in Eq. ^, integrate over 
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Figure 4: (Color online) (a) The effective action S{(j>) = 
S[r,(j)]/Vg as a function of for A = u/g = 100 and vari- 
ous fo- The different curves correspond to fo = 197, 195.94 
(fo'"), 195, 194, 192.9 (fo), 192, and 191 (from top to bot- 
tom), (b) The (red) solid curve shows the order parameter 
as a function Afo — fo ~ Vq. The (green) dashed curve 
shows the expectation value (j) of the metastable phase for 
ro < fo^. (c) The inverse susceptibility of the singlet mode 
1/Xs (X f — 2(j} as a function of Afo. 



(Pq, and absorb the factor T into the couplings. We ob- 
tain 



3V 



+ 



Stt 



6 It 
{r - 20) log 



r-20 



3r 

+ 2(r + (/))log 



(8) 



A2 



where A is the upper momentuni cutoff, and tq, which 
is proportional to the temperature, is the control param- 
eter. The renormalized r is proportional to the inverse 
magnetic correlation length £^~^. Differentiating Eq. ([8|) 
with respect to r and (j) yields the coupled saddle-point 
equations which determine their values at a given ro: 



ro 



2 log 



A6 



(r-2(/))(r 



(9) 



where fo = ro/g + {3u/2g)\og{A'^/g) ex (T - Tat), and 
T/v is the rescaled mean-field T^. The ratio A = u/g 
is large in our model, of order W/Tn, where W is the 
bandwidth [l3|. The analysis of Eq. ((TTl) for A » 1 
shows that the first non-zero solution appears at a par- 
ticular temperature when f™ w | A log 3 and at a finite 

(f) w 2.15 + 14.2/A. This obviously indicates that the 
mean- field Z^ transition is first-order. The actual transi- 
tion temperature is smaller than fp™ because at fp" the 
effective action only develops a local mininmm at nonzero 
(j), but this may not be a global mininmm. To find when 
the actual transition occurs, we solve Eq. (jH]) for r{(j)) nu- 
merically, substitute the result into ^ and obtain the 
effective action S{4>) for which fg is a parameter and 
Eq. pn]) is the saddle-point solution. The behavior of 
S'(0) for various fo is shown in Fig. Ul^a). At sufficiently 
large fo, it increases monotonically with and its only 
minimum is at = 0, implying that Z^ is unbroken. At 
fo = rj", the function S{(f)) develops an inflection point, 
which at smaller fo splits into a maximum and a mini- 
mum. At some fo = Tq the value of S'(0) at this minimum 
becomes equal to S{0), and for fo < Tq , the global mini- 
mum of the free energy jumps to a finite 0^0. Once this 
happens, the system spontaneously chooses one out of 
four states with ±0, and the Z^ symmetry breaks down. 
We plot (f) versus fo in Fig. [Hb) . 

To find how much the Z4 transition temperature Tz^^ 
actually differs from the SDW transition temperature 
T/v, we computed the spin susceptibility xil) within 
RPA, explicitly related fo to (T — Tn), and expressed A 
in terms of the ratio of Tjy and the fermionic bandwidth 
W. Collecting all factors we find 



Tz 



Tm + a 



-'- N 



1 



W \ogW/T, 



N 



(12) 



= ff log 



r-20 



(10) 



where the rescaled parameters are u = 3u/47r and g — 
g/4:Tr. For = (i.e., when Z4 is unbroken), Eq. ^ 
takes the familiar form r + 2>u/ (21ogA2/r) ~ tq whose 
solution is positive for arbitrary ro, i.e. 0(3) symmetry 
is preserved at any non-zero T. This is an obvious con- 
sequence of the Mermin- Wagner theorem. However, the 
discrete Z4 symmetry can be broken at a finite T. 

Assuming that is non-zero and eliminating r from 
Eqs. PT?)l and ^, we obtain the self-consistent equation 
for the rescaled order parameter = 0/g: 



2e^ + l 
e-^- 1 



+ A 




ro, (11) 



where a = 0(1), and T/v is the "mean field" Neel temper- 
ature, which does not take into account the suppression 
of SDW order by thermal fluctuations. The actual T)v 
tends to zero in 2D, but Tz^ remains finite. 

To analyze how the broken Z4 symmetry affects SDW 
correlations, we compute the eigenvalues of the spin sus- 
ceptibility matrix Af in ([7]). The two eigenvalues corre- 
spond to a singlet and a doublet mode Xs = l/('' — 20) 
and Xd = l/l?" + </')■ If either r — 20 or r -|- jumped to 
a negative value at the Z4 transition, then the breaking 
of Z4 would induce a simultaneous breaking of the 0(3) 
symmetry. However, it follows from pO)l that both Xs 
and Xd remain finite when jumps to a nonzero value, 
i.e. breaking the Z4 symmetry does not induce SDW 
order immediately (see Fig. [He)). 



C. Beyond mean-field: 4-state Potts model 



The effective action 
and large A as: 



can be expanded for small 



V 



(ro 



ron4>'-^J' + ^J' 



12 



16 



(13) 



This action has the same form as that of the 4-state Potts 
model J21|, implying that both transitions belong to the 
same universality class. We can use this analogy to go 
beyond the saddle-point solution and understand how the 
Z4 transition is affected by fluctuations of <j> fields. The 
4-state Potts model in 2D does exhibit a transition, i.e. 
the preemptive Z4 ordering is not destroyed by fluctua- 
tions ^22i] . Interestingly, however, fluctuations transform 
the first-order transition into a second-order transition, 
although with a rather small critical exponent (3 = 1/12 
for ~ (Tc - T)f^ (Ref. [Hi). A small /3 implies that 
the order parameter sharply increases below the critical 
temperature, and in practice this behavior is almost in- 
distinguishable from that in the first-order transition. 

Notice that the cubic term in the action, which comes 
from the product 0i(/)203, ensures that the Z4 symmetry 
breaking belongs to the universality class of the 4-state 
Potts model, and not of the 4-state clock model. This 
distinction is important, as they have different critical 
behaviors in two dimensions. While the 4-state Potts 
model has a /3 = 1/12 exponent, as discussed above, the 
4-state clock model transition belongs to the same uni- 
versality class of the Ising model, with /3 — 1/8. 



D. Experimental manifestations 

As spin rotational symmetry is preserved in the pre- 
emptive Z4 phase, no magnetic Bragg peaks are to be ob- 
served in neutron scattering experiments. On the other 
hand, since the charge density p{r) and the Casimir oper- 
ator S^(r) have the same symmetry, a spatial modulation 
of the latter induces a modulation in the charge density. 
Given the 2D character of graphene, such a super-lattice 
structure can be directly probed by STM. The additional 
Bragg peaks due to the quadrupled unit cell should also 
be detectable by scattering measurements. Local probes 
such as NMR can measure the different on-site fluctu- 
ating magnetic moments of the Z4 phase, since the size 
of the local moment controls the linewidth of the NMR 
signal. We thus expect to see two different linewidths 
coming from the 3A and the A sites. 



The order parameter can also be inferred by mea- 
suring the static magnetic susceptibility x ^-t any of the 
three nesting vectors. In the absence of 0(3) breaking, 
we have x(^o) = (2xd + Xs)/3- Once the order param- 
eter (j) jumps to a finite value below the transition, so 
does the susceptibility x(ro) = r^^ + (^^f^'' + ■ ■ ■ , where 
f is the value of r at = 0. This provides a direct 
method for detecting the order parameter 0. The jump 
of the static susceptibility (i.e. of the spin correlation 
length) also affects the electronic spectrum. For larger 
correlation length the system develops precursors to the 
SDW order, which give rise to a pseudogap in the elec- 
tronic spectral function. This pseudogap can be probed 
by photoemission experiments [6|. 



IV. CONCLUSION 

We discussed in this work the intriguing possibility of 
an emergent paramagnetic phase with spontaneously bro- 
ken translational symmetry for properly doped fermions 
on triangular and hexagonal lattices. This unique state 
emerges from a preemptive phase transition which breaks 
only a discrete translational Z4 lattice symmetry but pre- 
serves 0(3) spin- rotational invariance. We demonstrated 
that this phase exists in 2D systems and by continuity 
should exist in anisotropic 3D systems. We argued that 
such a phase should be observed in STM, NMR, neutron 
scattering, and photoemission experiments. 
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Appendix A: Hubbard-Stratonovich transformation 

In this Appendix we present the details of the 
Hubbard-Stratonovich transformation for the preemp- 
tive phase. We first introduce six bosonic fields ^/j, (^1, 
C2, and (/>i {i = 1,2,3), each corresponding to one of 
the fourth-order terms in Eq. [S] of the main text. Ex- 
plicitly, the interaction terms in the partition function 
Z = J VA.i exp(— 5'[Ai]) can be rewritten as 



exp 
exp 



u 

72 


[{Aj + Al + Alf 

J X 


= / Vij) exp / 


g-^(A? + A^ + A 
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exp 
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= JVC2 exp / 
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exp 


[f/'-'-^'i 


= /"p0fccxp /" 


-0 + 0fc(A,.A,) 
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(Al) 
(A2) 
(A3) 

(A4) 



where {ijk) in the last equation are cychc perturmations 
of (123). The new action becomes 

5[A„V,C0] = 5^y"A'„[V',C,0](A,-A,) 

A V 2g 2i; 2u; ^ ' 

where C, — (Ci,C2) and 4> — {4'It4'2t4'3)^ and the matrix 
X is 



<%- = ! -03 X^1+2U2-C 

-02 -01 



with x„ ^ = ro+V'+q^ = ^+9^1 and the three miit vectors 
are Ui^2 = (1/2, ±V3/2), and U3 = (—1,0). Integrating 
out the Ai fields yields an effective action 



5[V',C,0] = -y log(deti[^sC,0]) 




2g 



2v 2u 



(A7) 



As discussed in the main text, the mean-field solution of 
the potential preemptive phase is given by the saddle- 
point solution of coupled equations: dS/dip — dS/dC, = 
dS/dcj) = 0. In particular, we consider the two equations 
involving the doublet Q: 



Ci 

C2 



3w 

y 
y 



6CI - 6C1 (x,-^ + Ci) - (0? + 



201) 



6C2 (2Ci 



dctAf 



detA- 



(A8) 



It can be easily checked that the mean-field configura- 
tions with Ci = C2 = and |0i| = |02| = |03| = are 
solutions of the above two equations, indicating that the 
Z4 phase solutions discussed in the main text satisfy the 
saddle-point equations of the effective action (|A7p . 
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